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Abstract. We present a new fast method for simulating pencil-beam type light-cones, using the MLAPM-code 
(Multi Level Adaptive Particle Mesh) with light-cone additions. We show that by a careful choice of the light-cone 
orientation, it is possible to avoid extra periodicities in the light-cone. As an example, we apply the method to 
simulate a 6 Gpc deep light-cone, create the dark matter halo catalogue for the light-cone and study the evolution 
of haloes from z = 6 up to the present time. We determine the spatial density of the haloes, their large-scale 
correlation function, and study the evolution of the mass function. We find a surprisingly simple relation for the 
dependence of halo maximum mass on redshift, and apply it to derive redshift limits for bright quasars. 
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1. Introduction 

The recent (completed and on-going) galaxy redshift sur- 
veys (e.g. LCRS, 2dFGRS, 6dF, SDSS) of our local uni- 
verse have enormously increased our knowledge of galax- 
ies and of their large scale distribution. These large-scale 
structure data sets provide also a powerful tool for break- 
ing many of the parameter degeneracies associated with 
the CMB data (Spergel et al. I2003[) . Clusters of galaxies, 
the largest virialised systems known, have a vital role in 
understanding cosmological structure formation. In par- 
ticular, following the evolution of clustering with redshift, 
we can put direct constraints on models for the evolution 
of density perturbations (Munshi ct al. 2004). 

However, it is impossible to study cosmological evo- 
lution, using only such local surveys. Thanks to recent 
developments in instrument technology there are several 
new surveys designed to probe at increasingly higher red- 
shifts, in order to study the early evolution of galaxies and 
their systems. We shall list a few of them. 

The DEEP2 redshift survey is planned to study the 
evolution of properties of galaxies and the evolution of 
the clustering of galaxies from z ~ 1.5 to z = 0, using the 
DEIMOS spectrograph on the 10 m Keck II telescope. The 
results will be compared with those of the local surveys, 
e.g., as the LCRS (Coil et al. 12002) ■ The ALHAMBRA 
photometric survey (Moles et al. I2005|l has similar goals. 
It will cover 4 square degrees of sky, and will find galaxies 
up to redshifts z ps 5. 
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The VIRMOS galaxy redshift survey, using VIMOS 
on the VLT telescope, is also designed to study the for- 
mation of galaxies and large-scale structure over the red- 
shift range (0 < z < 5), covering sixteen square degrees 
of the sky in four separate fields (Virmos Consortium, 
http: / /www. oamp.fr/virmos/ ). 

GOODS (the Great Observatories Origins Deep 
Survey) aims to unite extremely deep multi-wavelength 
observations (up to z ~ 6) both from ground-based 
(VLT, Keck, Gemini, NOAO, Subaru) and space tele- 
scopes: Hubble (HST Treasure program), SIRTF (Legacy 
Program), and Chandra. Observations cover two fields 
(10' x 16'), centered on the Hubble Deep Field North and 
the Chandra Deep Field South. The primary goal of the 
GOODS program is to provide observational data for trac- 
ing the mass assembly of galaxies throughout most of cos- 
mic history (Dickinson & Giavalisco 2003). 

The Galaxy Evolution Explorer (GALEX) satellite was 
launched in 2003 (Martin et al. 2005 ) . In combination with 
ground based optical observations this satellite allows to 
study the star formation rate, the galaxy luminosity func- 
tion and other parameters of galaxy evolution over the 
redshift range < z < 5, and even beyond this limit (see 
http://www.galex.caltech.edu). 

Several recent studies have revealed clustering of Ly 
alpha objects, quasars and radio galaxies at very high red- 
shifts (Ouchi et al. I2UU41 Venemans et al. I20U1 Malhotra 
et al. IMiH Wang et al. l!?UU5l a nd Stiavelli et al. Jfflty . 
For example, Ouchi et al. ( 2004) reported about the dis- 
covery of primeval dense structures at redshift 6 contain- 
ing Lyman alpha emitters (LAEs). Ouchi et al. estimated 
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that masses of these structures - progenitors of present-day 
large scale structures are of order of 10 12 Mq - 10 13 Mq 
and dimensions are less than 10 Yr 1 Mpc. 

By these surveys, observations are reaching the scales 
where the evolution of galaxies along the light-cone plays 
an essential role. Comparison of deep surveys with our 
theoretical understanding of structure formation and evo- 
lution has become even more topical and thus simulations, 
or mock catalogs, have become an important tool for the 
design of observational projects and for later data anal- 
ysis. A properly constructed simulation can allow us to 
compare the observations directly with modern models of 
structure formation, can serve to test the data process- 
ing algorithms for biases, and to quantify the impact of 
numerous observational selection effects (Yan et al. 2004). 

We believe that light-cone simulations mimic better 
the observational datasets than conventional simulations 
at z = 0, and at different "snapshots" at fixed redshifts. 
The light-cone algorithm stores 'light-cone particles' on- 
the-fly, allowing us to follow continuously the evolution 
of structure. Although not difficult in principle, light-cone 
simulations tax heavily computer resources - the volumes 
are too large to model them with necessary mass reso- 
lution. The best-known light-cone, "the Hubble simula- 
tion" (Colberg et al. 12000(1 . was run on specialised parallel 
shared-memory computers, and remains the only one with 
publicly available data. 

The goal of this paper is to present a new method 
to perform light-cone simulations of deep pencil beams, 
typical for deep surveys. Our method is rather efficient 
and allows to perform and analyse the simulation on single 
workstations. We shall build a light-cone from z = to 
z = 6 in a "concordance" cosmological model, and shall 
follow the evolution of dark matter (DM) haloes with the 
cosmic epoch. Visualizations of DM-haloes in the light- 
cone model can be seen at Tartu Observatory web pages 
(http : //www . aai . ee/~maret/lc .html). 

When our project was almost finished we learned that 
a similar light-cone project MoMaF has been realized by 
Blaizot et al. 1(200511 . Below we shall compare their ap- 
proach with ours. 

2. The light-cone simulation 

2.1. Method 

We use the Multi Level Adaptive Particle Mesh code 
(MLAPM, Knebe et al. lMiTjl with light-cone output. The 
MLAPM code is adaptive, with sub-grids being adap- 
tively formed in regions where the density exceeds a spec- 
ified threshold. There are similar codes, as the ART code, 
Kravtsov et al. (J1997); it is essential that simulations are 
done with different codes, to cross-check the results of the 
complex dynamics of the cosmic structure formation. 

We base the light-cone on periodic replicas of our sim- 
ulation cube. This could, of course, introduce extra peri- 
odicities in the results, but we shall demonstrate that a 
clever choice of the light-cone parameters will avoid those, 



at least for pencil-beams. If we choose the pencil-beam 
along one of the coordinate axes, we will sample the same 
regions of the cube many times; but if the direction of 
the beam is oblique, and the beam narrow, we will use 
different regions of the cube along the beam. 

In order to collect the light-cone, we find first the 
co-moving light-cone radius and determine all the copies 
of the original simulation cube, which intersect with the 
spherical surface of that radius (cover it). For the sim- 
plest geometries (almost full-sky) that is enough, but for 
a pencil-beam patch geometry we have to use three addi- 
tional checks: 

1. find if the cube copy vertices lie inside the patch; 

2. find if the patch edges intersect the copy's faces; 

3. find if the patch planes intersect the copy's edges. 

If nothing this happens, the copy-cube does not in- 
clude the light-cone. This cube list may seem to be an 
unnecessary complication, but it will speed up light-cone 
processing for the most interesting deep and narrow light- 
cones. 

When the cube list is completed, we check for all par- 
ticles their positions in all the candidate cubes, and if the 
(copy) particle has crossed the light-cone between the pre- 
ceding z and the present z, we find the crossing redshift 
and coordinates by linear interpolation. Then we check 
for the patch geometry, if necessary, and write the (copy)- 
particle data. 

Detailed instructions on generation of the light-cone 
simulations can be found in the MLAPM user guide 
Hhttp: / /www.aip.de/People / AKnebe /MLAPM|) . 

The MoMaF project fBlaizot et al. 12005(1 takes a 
slightly different approach, saving first a large number of 
snapshots and using these later to build the light-cone. 
When assembling the light-cone, they randomly rotate the 
snapshots, in order to eliminate the influence of periodic- 
ities. As this destroys the continuity of the structure, we 
prefer not to scramble our cubes. In fact, MoMaF allows 
such an approach too, for pencil-beam surveys, although 
Blaizot et al. (2005J devote most of their attention to the 
scrambled case. 

2.2. Application 

We simulate a pencil-beam mock survey, similar to the 
contemporary observational surveys. We use a simulation 
with 256 3 dark matter particles in a 256 3 grid. We follow 
the evolution from z = 6 to present time. Our simulation 
covers 2° x 0.5° in the sky. 

We use a flat cosmological model with the parameters 
derived by the WMAP microwave background anisotropy 
experiment (Bennett et al. 12003(1 : the dark matter density 
fi m = 0.226, the baryonic density = 0.044, the vac- 
uum energy density (cosmological constant) Q\ = 0.73, 
the Hubble constant h — 0.71 (here and throughout this 
paper h is the present-day Hubble constant in units of 
100 km s _1 Mpc" 1 ) and the rms mass density fluctuation 
parameter erg = 0.84. 
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The initial data for our simulation was gener- 
ated in a cube of 200 hr 1 Mpc co-moving size, for 
z = 30. Thus, each particle has a mass of 3.57 • 
10 10 /i _1 Mq. The transfer function for our model was 
computed using the COSMIC code by E. Bertschinger 
(http : / / arcturus .mit . edu/cosmics/). 

Figure ^ shows the projected (dark-matter) density 
fields of our light-cone simulation from z — to z = 6 
(for the cosmological model chosen z = 6 corresponds 
to 5981.42 h~ 1 Mpc. The density was calculated, using 
an Epanechnikov kernel (h(x) ~ (1 — x 2 /r 2 ) of radii r 
1.0 hr l Mpc, 16.0 ft, -1 Mpc - i.e., for the characteristic 
scales of clusters and superclusters of galaxies. 

The full amount of light-cone data (positions and ve- 
locities of dark matter particles) can have enormous. The 
evident way to compress the data is to extract catalogues 
of dark matter haloes; this are the objects that can be 
associated with observed galaxies and clusters of galaxies. 

As the amount of data is large, we used the simplest al- 
gorithm (friends-of- friends, FOF) to identify dark matter 
haloes. The simulations reported in this paper were run 
in early 2004; now we would have used a new feature of 
MLAPM, its halo finder (see Gill et al. l2"UU4) that uses the 
basic adaptive density tree of the simulation, and outputs 
haloes at all required snapshots. It should not be difficult 
to adapt that to output light-cone halo catalogs on the fly; 
this remains to be done. Deep light-cones generate huge 
amounts of data, and the only practical way to handle this 
is to keep only the haloes. 

FOF uses a linking length b to collect particles in 
groups with spacing closer than b times the mean inter- 
particle spacing. The linking length is frequently deter- 
mined from the virialisation density p v (see Bryan and 
Norman 1998), obtained from the solution for the collapse 
of a spherical top-hat perturbation. The density A c depen- 
dence on background cosmology via the matter density 
parameter 

n M (z) = n (i + zf/E(z) 2 , 

where 

E 2 = fi„(l + zf + R (1 + z) 3 + n A . 

Bryan and Norman 1998. give an approximate formula for 
that: 

A c = Af/Pcrit = 1871-2 + 82x ~ 39a;2 ' 

where x — Qm(z) — 1 and p cr it is the density in the 
Einstein-deSitter model (0 A = 0, fijvf = 1). The actual 
mean density is Q(z)p cr j±. Thus the linking length that 
would select virialised haloes at a given redshift z can be 
written as b = (£Im ( Z )A C ) 5 • 

The assumption that objects populate only virialised 
haloes could be a bit extreme. We normalised the link- 
ing length by its present value, comparing the amplitude 
of the simulated halo mass function (at it's massive end) 
with an observational mass function of galaxy groups of 
the Las Campanas Loose Groups of Galaxies, hereafter 



LCLG, Heinamaki et al. (2003). In order to obtain this 
mass function, we ran the MLAPM code in the snapshot 
mode, with the same initial data we used for the light-cone 
simulation. 

Using this normalisation, we chose the present epoch 
(z = 0) linking length as b — 0.23 (in units of the mean 
particle separation) , which corresponds to the matter den- 
sity contrast 8n/n = 80. This is somewhat lower than the 
virial density, but not much, and delineates well galaxy 
groups, as seen in the LCLG Catalog (Tucker et al. 2001). 

Figure |3shows how the linking length we used changes 
with redshift. Due to the accelerated expansion at later 
epochs the haloes need higher density contrasts (and 
smaller linking lengths) for virialisation. 
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Fig. 2. The FOF linking length as a function of redshift, 
for the light-cone simulation. 

We extracted several populations of dark matter haloes 
from our light-cone: haloes with 2 or more particles (Mh > 
7.1 • 10 10 /i _1 Mq), haloes with more than 20 particles 
(Mh > 7.1 • 10 /i Mq), and a conservative sample of 
haloes with more than 100 dark matter particles (Mh > 
3.6 • 10 12 /i. -1 Mq), We cleaned all the group catalogues, 
using the virial equilibrium condition E r = Eki n / \ E pot \ < 
0.5 (Ep t is the potential energy and Ekin - the kinetic 
energy of a group) for groups to be included in our final 
group catalogue. 

The number of small haloes is considerably reduced by 
the virial condition, which is natural - many of such haloes 
are random encounters. From the ^ mm — 20 halo sample, 
6.5% of haloes are rejected, and only a few haloes from the 
rich halo sample do not satisfy the virial condition. This 
indicates that the linking length we have chosen does not 
pick up too many non-virialised objects. The final full halo 
sample includes about 600,000 haloes, the intermediate 
sample has 14016 haloes and the sample of rich haloes 
includes 4088 haloes. 
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Fig. 1. The density field of the 6 Gpc deep dark matter light-cone simulation smoothed with a r = 1.0 h 1 Mpc kernel 
(upper panels), and with a r = 16.0 hr x Mpc kernel (lower panels). 



3. Checking for periodicities 

Assembling a light-cone from periodic replicas of a simula- 
tion cube will lead to several statistical biases. These are 
thoroughly discussed in Blaizot et al. (200SJ- The main 



bias comes from scrambling the snapshots; as we do not 
use scrambling, we bypass this bias. 

Another bias is caused by the finite volume of the sim- 
ulation cube and can be corrected for; Blaizot et al. (2005) 



P. Heinamaki et al.: LCS 



5 




0.9 - * 



0.8 I ' ' ' ' ' 1 

50 100 150 200 250 300 

r (Mpc/h) 

Fig. 3. The large-scale correlation function of DM haloes 
in the light-cone. The solid line with error-bars shows 
the correlation function for the whole sample, clotted bars 
show the correlation function and its 1-a deviation for the 
redshift interval z e [0, 1]. 

show how to do that. We can also, in principle, to modify 
the light-cone by linear large-scale modes; this would al- 
low us to model better correlation functions and number 
counts. 

There is one more effect to worry about - as our light- 
cone is composed of periodic replicas of a much smaller 
simulation cube, it could have periodicities at scales of 
about the cube size. The easiest way to check for that is to 
look at the two-point correlation function of the light-cone. 
The total number of light-cone mass points is very large; 
in order to calculate the correlation function we chose the 
full halo sample, and diluted it randomly, by a factor of 
six, to keep the 0(N 2 ) calculation time within reasonable 
limits. The number of haloes used was about 10,000. We 
are worried about large-scale correlation, thus we chose 
the Landy-Szalay estimator, which is one of the best es- 
timator for large scales (see, e.g., Martinez & Saar l2002|l . 
The rms error of the estimator was taken to be twice the 
Poisson error. We show this correlation function in Fig. El 

As we see, no periodicity at scales around 200 h^ 1 Mpc 
can be seen. The dip in the correlation function at ~50 
h~ 1 Mpc and the maximum around ^120 ft, -1 Mpc are 
real, and are seen both in simulations and in the observed 
large-scale correlations (e.g., for Abell clusters, Einasto 
et al. 119971 see also Einasto et al. 119941 and Tago et al. 
2002). A slight indication of periodicity could be seen in 
the closest redshift interval (z = ... 1, shown with dot- 
ted error-bars). However, even there £(r) = for distances 
r > 30/i _1 Mpc, within 2a limits (the volume and the 
number of haloes are small, and the estimate is rather un- 
certain). For other redshift intervals the correlation func- 
tions practically coincide with the mean for the whole sam- 
ple. 

This can be explained by the good choice of the light- 
cone geometry. An oblique direction of the light-cone with 
respect the symmetry axes of the cube will cause the light- 
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Fig. 4. The volume multiplicity histograms of the light- 
cone. 

cone to sample different regions of the simulation cube in 
its different replicas. We can quantify this by counting the 
number of times each cube volume element is included in 
the light-cone. In order to estimate that, we choose a sim- 
ple algorithm - we populate the light-cone with a Poisson 
point process, fold it back onto the original simulation 
cube, and find the resulting one-point density distribution 
in this cube. When properly normalised, this distribution 
approximates the volume multiplicity distribution. 

The volume multiplicity distribution for our light-cone 
is shown in Fig.0J The two different histograms show the 
results for two different oversampling factors (the ratio of 
the density of the Poisson process to the mean light-cone 
density) - the higher this factor, the better the estimate, 
and the larger the computer time. 

We see that most of the volume of the simulation cube 
is used only twice in our light-cone. The worst choice of the 
light-cone direction, along an axis of the cube, would give 
a factor close to 30 (the depth of the light-cone divided by 
the cube size). We can use the mean of the histogram or 
its rms error as the figure-of-merit M for the volume use; 
maybe the best choice would be to use a a combination of 
these: 

M 2 ^a^+fh 2 , 

where m is the volume multiplicity, and fh its mean value. 
Thus, the figure-of-merit for an ideal case (e.g., the Hubble 
volume 7r/8 sterradian light-cone, Colberg et al. I2000f) 
would be 1; for our light-cone it is 2.9 (but could be over 
30). Another relevant characteristic could be the distri- 
bution of distances between point replicas. A rough es- 
timate of the mean replica distance for our light-cone is 
6000/3=2000 h^ 1 Mpc; this is why we do not see period- 
icities in the correlation function. 

The light-cone figure-of-merit can be found before sim- 
ulations, and can be used to select the best orientation 
parameters for a given light-cone geometry (its size and 
the size of the simulation cube) . Its value will also charac- 
terise the possible periodicities in the light-cone. We plan 
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Fig. 6. Spatial density of the haloes as a function of red- 
shift for five different mass limits. 



Fig. 5. The density field of the simulation. Halos are 
marked with black circles. The figure shows the most dis- 
tant part of the light-cone, from 4 to 6 Gpc. The up- 
per panel shows rich haloes (M > 3.6 • 10 12 /i _1 Mq), and 
the lower panel shows also the intermediate-mass haloes 
(M > 7.1 • lO n /i- 1 M ). 

to include the figure-of-merit code in the tools section of 
the MLAPM code package. 

4. Evolution: redshift dependence of the dark 
matter haloes 

4.1. Spatial density of haloes 

Light-cones as ours can be used for many purposes. As this 
paper is mainly meant to present the light-cone construc- 
tion method, we shall give only a few examples below. 

First, as our light-cone is pretty deep, we can follow 
the formation of first massive haloes. Smaller haloes form 
earlier yet, and their study would require a deeper light- 
cone. 

Figure |3] illustrates the appearance of the first dark 
matter haloes. In panel (a), DM haloes with masses M > 
3.6 • 1O 12 M0 are shown. Halos, placed on the density field 
map (smoothed with 1 hr x Mpc kernel), are marked with 
black circles. As we see, the very first massive haloes ap- 
pear at about D=5600 h^ 1 Mpc (z=4.97). Panel (b) shows 
shows haloes of mass M > 7.1 ■ 10 11 Mq. In this case the 
very first haloes appear already at D = 5900 h^ 1 Mpc 
(z = 5.62). Smaller haloes exist in our light-cone from its 
start, z — 6 already. 

The spatial densities phaio of all our DM haloes as a 
function of redshift (the number of haloes per co-moving 
volume at redshift bins of width 0.2) for different halo 
mass limits are shown in Figure Each line indicates 
different halo mass limits, from 10 10 Mq (at the top) to 
10 14 Mq (at the bottom of the figure). As can be ex- 
pected, rich haloes appear later than poor ones. Haloes 



with masses smaller than 10 12 Mq exist in our light-cone 
from the beginning. Overall density evolution of haloes 
of the two smaller mass ranges is the same; in the red- 
shift range z = 5.7. . .4 the spatial density of haloes in- 
creases about an order of magnitude. More rapid increase 
occurs in the redshift interval z — 3.7. . .3. These transi- 
tions are also seen in the density behaviour of haloes of 
masses M > 1O 13 M . 

The halo density reaches its maximum at ~ z = 0.9, 
after that it decreases. It could be caused by merging of 
smaller haloes, but could also be a statistical fluctuation, 
due to the small volume of our pencil beam at redshifts 
less than one. 

The first haloes of mass M > 10 14 Mq appear at z ~ 
3. The large variation in their density evolution is due 
to rareness of such objects and to a small volume of the 
sample. There is also two overall slopes seen in mass scales 
less then 10 14 Mq. We also notice that density evolution 
on almost all mass scales (excepting the the largest mass 
interval) between z=5.5 and z=3 is more rapid than the 
later density evolution, between z=3 to z=l. The reason 
for that could be gradual merging of smaller haloes to 
form bigger systems. 

4.2. Halo mass history 

Figure shows the cumulative mass function (MF) for 
dark matter haloes, defined as the number density of 
haloes above a given mass M, n(> M). The solid line 
shows the mass function for the intermediate mass light- 
cone haloes (n > 20). The dotted line shows the MF for a 
similar halo sample, obtained in our MLAPM simulation, 
for the final snapshot (a 200 Mpc box at z = 0). The 
mass functions differ by an order of magnitude at 10 11 Mq, 
due to the smaller overall halo density in the deep light- 
cone. Another difference can be seen at the high mass end 
of the MF-s - the z = MF extends to larger masses. This 
is due to the narrowness of the pencil beam, that occupies 
only about 10~ 4 of the simulation box at z = 0. Statistics 
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Fig. 7. The halo mass function for the light-cone (inter- 
mediate sample, n > 20) (solid line) and for a conventional 
simulation for the present time (dashed). The dotted line 
shows the light-cone halo mass function for all haloes. 

of such surveys should be taken with caution, at least for 
redshifts less than z sw 1. Without this cosmic variance 
effect, the light-cone MF should follow the snapshot MF 
throughout the whole mass range. The dotted line shows 
the light-cone mass function for all halo masses; this grows 
steeply, showing that small haloes dominate in number, es- 
pecially at early epochs, when massive haloes did not exist 
yet. 

The light-cone can be also used to study the evolution 
of the halo mass function. We plot this function for dif- 
ferent redshift intervals in Fig. El The overall shape of the 
mass functions between z — to z — 1 and z = 1 toz = 2 
is quite similar to the mass function of the whole sam- 
ple, for the mass interval M > 1O 12 M (see Fig.0). The 
mass functions steepen with redshift, showing the domi- 
nance of small haloes at large redshifts. The characteristic 
halo masses decrease with redshift, with the maximum 
halo mass at the largest redshift interval z = 5 ... 6 being 
1.5- 1O 12 M . 

The most easier objects to observe at any epoch are 
the most luminous (most massive) ones. In Fig.[§]we plot- 
ted the maximum halo mass for a given redshift. This 
figure shows that the maximum masses increase exponen- 
tially with decreasing redshift, down to the redshift ~ 0.9, 
where the cosmic variance limit appears. This interesting 
trend, shown in the figure, can be approximated by an 
intriguingly simple formula 

M max (z) « 10 15 M Q 10-°- 5z « 10 15 M Q e- z . (1) 
This relation is shown by a dashed line in the figure. 

5. Early haloes and black holes 

As an example of an application of light-cone models, we 
discuss an interesting possibility that masses of supermas- 
sive black holes (BH) may correlate with masses of dark 
matter (DM) haloes. 



Fig. 8. Halo mass functions for different redshift intervals. 



Using a sample of 37 local galaxies Ferrarese (2003 ob- 
tained a relation between the masses of BH and circular 
velocities of galaxy disks, similar to the well-known rela- 
tion between the' BH mass and the bulge velocity disper- 
sion (e.g. Magorrian et al. lTMSl M BH /M bu i ge = 3 • 10~ 3 ). 
Using Bullock et al.'s (200 1 • simulations (their result 
M D m = 2.7- 10 12 (?w/200 km s _1 ) 3 M agrees well with 
our results) together with Ferrarese's observational re- 
sults, the best fit for the ratio between Mbh and Mom 
can be written as 

M BH /10 8 M Q ~ 0.10(M DA/ /10 12 M Q ) 1 - 65 (2) 

(see figure 5 in Ferrarese's paper). Ferrarese pointed out 
(see the references therein) that the critical DM mass for 
the BH formation could be about 5 • 1O 11 M , while less 
massive haloes might not offer favorable conditions for BH 
formation. According to the previous equation, this limits 
the BH masses by ~ 1O 6 M in such DM haloes. Ferrarese 
concluded also, that according to the results by Hachnelt 
et al. I|1998f) only haloes with masses M > 1O 12 M may 
host a black hole with Mbh > 1O 6 M . In our light- 
cone the spatial density of such host haloes is less than 
2 ■ 10" 7 Mpc~ 3 /i 3 at redshifts larger than z = 5.5. 

If we assume that Ferrarese's Msn/MDMhaio relation 
also holds at high redshifts, and combine it with our max- 
imum halo mass relation l[T]l. we obtain a simple relation 
for the BH masses likely to reside in the most massive 
dark matter haloes as a function of redshift: 



M 



BH 



10 12 M 







■(*) = 10 



-0.83z 



(3) 



McLure & Dunlop (2004) concluded, based on virial 
BH mass estimates from the SDSS first release and for 
12698 quasars, that quasars' SBH (super massive black 
hole) masses lie between ~ 1O 7 M and ~ 1O 9 M . If we 
substitute these limits in equation JSJ, we find that the 
most prominent quasars exist between the redshifts ~ 6 
to ~ 2. Haloes massive enough to form SBH and hence 
quasars, either did not exist before the redshift z rs 6, 
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Fig. 9. Masses of the most massive haloes for different 
redshifts 



or they were very rare, with the number density < 
W- 7 Mpc- 3 h 3 . 

6. Conclusions and perspectives 

We have presented a new method for simulating pencil- 
beam type light-cones, using periodic replicas of a base 
A^-body cube. Such light-cones are typical for extremely 
deep optical surveys, either underway or in the planning 
stage. We have shown that by a careful choice of the light- 
cone parameters, it is possible to avoid extra periodicities 
in the light-cone. 

We have simulated a deep (up to z = 6) light-cone, 
generated the dark matter halo catalogue, and studied its 
properties. We find that early light-cone is dominated by 
small haloes, and the maximum halo mass can be clearly 
traced throughout all the epochs. 

We find a simple approximation for the dependence of 
maximum dark matter halo mass on redshift, and use it 
to explain the redshift limits of the quasar distribution. 

In summary, the algorithm we use to simulate light- 
cones is lightweight and fast. The code is in the public 
domain, included in the MLAPM package. Nevertheless, 
the main bottleneck of any light-cone model is the vast 
amount of data, proportional to D 3 (D is the depth of the 
light-cone). The only way to overcome this is to discard 
most of the hard-obtained simulation data, the positions 
and velocities of dark matter particles along the light-cone, 
and to store only the data on dark matter haloes, creating 
and analysing them on-the-fly. This could happen soon, 
as our base MLAPM code already includes the MHF halo 
finder, which outputs halo catalogues for fixed-time snap- 
shots. Such halo catalogs can be directly used to predict 
the SZ effect from early galaxy clusters for the Planck mis- 
sion. Moreover, as gravitational lensing directly traces the 
total matter density in the universe, our light-cone cata- 
logs are good tools to study the weak lensing of the light 



emitted by distant quasars, and of the CMB. These are 
the lines of our present work. 

Another important problem in simulating light-cones 
to compare with observational data is populating dark 
matter haloes with galaxies, and assigning the galaxies 
the features an observer sees. In fact, the MoMaF project 
(Blaizot et al. I2005|) already gives galaxy populations in 
light-cones, using semi-analytic recipes. We plan to in- 
clude such recipes in our code, also. 
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